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Abstract. We have detected gamma-rays from RXJ1713— 39 
in the energy range between 400 GeV and 5 TeV using a new 
kind of analysis: likelihood analysis. The statistical signifi- 
cance of the measurement is greater than 8 a. The details of 
this analysis method are presented. 



celeration by a shock wave produced by a supernova explo- 
sion and inverse Compton scattering with micro-wave back- 
ground radiation (Pohl , 1996; Mastichiadis, 1996a; Mas- 
tichiadis and de Jager , 1996b; Yoshida and Yanagita, 1997; 
Naito et al., 1999) can explain cosmic ray acceleration very 
well. 



1 Introduction 

The CANGAROO experiment is the Imaging Atmospheric 
Cherenkov Telescope located in Woomera, South Australia. 
It started with a 3.8-m telescope (Kifune et al., 1995). We 
now have a 10-m reflector (Kawachi et al., 2001) and are 
going to build a stereo-scopic system (Mori et al., 2001). 

Supernova remnants are one of the hot topics related to ori- 
gin of cosmic rays (Yanagita, 1999; Ellison et al., 2000). The 
CANGAROO experiment detected gamma-rays from SN1006 
(Tanimori et al., 1998a, 2001; Hara et al., 2001). Particle ac- 
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RXJ1713-39 was found in the ROAST aU sky survey (Pf- 
effermann and Aschenbach , 1996), and was found to have a 
shell structure. Hard X-ray emission was observed by ASCA 
(Koyamaet al., 1997). An association with a molecular cloud 
was found (Slane et al., 1999). The CANGAROO collabo- 
ration found an evidence for TeV gamma-ray emission from 
the northwest-rim (Muraishi et al., 2000). The preliminary 
result of CANGAROO-II (7m) observation in 1999 shows 
indication of a gamma-ray signal. In this report, an analysis 
of RXJ1713-39 with CANGAROO-II (10m) is presented. 
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Fig. 1. Shape parameter distributions: (a) distance, (b) length, 
and (c) width. The blank histograms were obtained from the OFF- 
source run. The hatched histograms are Monte-Carlo gamma-ray 
events. The cut positions are indicated by arrows. 



2 Analysis 

2.1 Data sample 
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The observation was carried out during two periods: 23-26 
July and 19-27 August, 2000. The pointing direction was as 
same as that of CANGAROO-I observation (Muraishi et al., 
2000), i.e., the NW-rim (RA, 5)={nHl'^hQ\7 , -39°31'52".4( J2000)), 
where the X-ray flux is maximum. The total observation pe- 
riods were 1419 and 1397 min., for ON- and OFF-source 
runs, respectively. We restricted the data to that taken at ele- 
vation angles greater than 60 degrees and without cloud, dew 
etc., by looking carefully at the proton shower rate. The re- 
sulting good quality data corresponded to 649 and 642 min. 
for ON- and OFF-source runs, respectively. After removing 
bad (hot) PMT and bright star hits, we applied a pulse height 
cut 3.3 photoelectrons) and a timing cut (±40nsec). Events 
with at least one cluster of five-adjacent triggered PMTs were 
then analyzed. 
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2.2 Conventional Imaging Cut 

First we rejected energetic multi-cluster events as follows. 
We calculated the energy-ratio {Edher / Emax) between the 
most energetic cluster (Emax) and the others (Eother)- Be- 
cause the gamma-ray events most likely have a single cluster, 
we rejected events with this ratio greater than 25%. Then, the 
event shape parameters (Hillas, 1985; Weekes et al., 1989) 
of distance, length, width, asymmetry and a were calcu- 
lated for the most energetic clusters (Fig 1). The blank his- 
tograms are obtained from the OFF-source run. The hatched 
histograms are Monte-Carlo gamma-ray events. The cut po- 
sitions are indicated by the arrows. Here, we refer to this as 
a square cut analysis. The cut dependences of the a dis- 
tributions are shown in Fig 2: Fig 2-(a) was obtained with- 
out any shape cuts, and (b) was obtained by a distance cut. 
By adding length and width cuts, we obtained Fig 2-(c). 
The data points with statistical error bars were obtained by 
ON-source runs, and the hatched histogram by OFF-source 
runs, respectively. Better signal-to-noise ratios (S/N-ratios) 
were obtained by tighter cuts. The normalization of OFF 



Fig. 2. Image orientation angle (a) distributions for ON-source data 
(solid squares) and OFF (hatched histogram), for a "square cut" 
analysis; (a) no cut, (b) distance cut, and (c) (distance, length, 
and width) cut. 



events was done using the number of events with a > 25 
degrees. This agreed with the time-interval ratios between 
the ON- and OFF-source runs. The obtained signal levels 
were: (a) 1510 ± 234(6.4cr), (b) 1129 ± 169(6.7cr), and (c) 
931 ± 127(7.30-), respectively. 

2.3 Likelihood 

It is well known that there are energy dependences in the 

standard shape parameters (Hillas et al., 1998). One way to 
minimize these effects is a likelihood analysis (Enomoto et 
al., 2001). Before doing this, we have observed a strong 
correlation between the distance and the other shape pa- 
rameters. This is considered to be due to the edge effect of 
the focal plane detector (i.e., camera). We, therefore, need 
a distance cut before starting a likelihood analysis. Be- 
cause the distance also has a small energy dependence, we 
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Fig. 3. Likelihood distributions for Monte-Carlo gamma rays (the 
hatched region) and OIT-source events (the blank histogram). 



Fig. 4. Image orientation angle (a) distributions for ON-source run 
(data points with error bars) and OFF (hatched histogram) for a like- 
lihood cut analysis. 



corrected it by a linear function of energy. After applying 

an energy-dependent distance cut, we calculated the likeli- allowing us to conclude that there are no event- specific bi- 

hood using the above shape parameters (length, width, and ases. The S/N-ratio was greatly improved by this analysis. 

asymmetry) for both gamma-ray and proton assumptions. The energy threshold for this analysis was estimated to be 

In order to estimate the gamma-ray probability-density-function 400 GeV with the present CANGAROO-II system. 
(PDF), we used Monte-Carlo gamma-ray events; for pro- The a distribution is wider than that for a point source 

tons, the OFF-source events were used. In the Monte-Carlo assumption (typically within 15°), consistent with the previ- 

simulation, gamma-rays with a Crab-Uke spectrum (E~'^-^) ous observation (Muraishi et al., 2000). Roughly speaking, 

were generated. In practice, we made two-dimensional his- the observed flux agrees with the previous observation (Mu- 

tograms, for example, log{ADC) vs length. We assumed raishi et al., 2000). A more precise flux calculation is now 

tiiat ADC is proportional to the energy of the showers. The being carried out, based on the following systematic error 

total number of events was normalized to unity and two- studies, 
dimensional (2D) PDFs were obtained. We defined the likelihood- 
ratio {L) as 

3 Possible Systematics 

L = Probin) / {Prob{-i) + Prob{p)) 

Here, we discuss the possible systematic errors in the accep- 
, where, Prob means product of 2D-PDFs of each shape pa- tance calculation using this likelihood method. The accep- 
rameter for the gamma-ray or proton assumption. Here, a tance was calculated using the Monte-Carlo method (Enomoto 
one-to-one contamination of gamma-rays to protons was as- et al, 2001). Gamma-rays with a Crab-like spectrum (E^'^-'^) 
sumed (this can be modified in the future), in spite of the fact were generated. The measured values of the detector param- 
that protons are dominant in cosmic-rays. The distributions eters, such as the point spread function (PSF) of the mirrors, 
of L are shown in Fig 3. The hatched histogram was L for the reflectivity, etc. we used. In order to check, we analyzed 
Monte-Carlo gamma-rays and the blank OFF-source events. Crab data (observed in December 2000 with the 10-m tele- 
The final event samples were obtained by cutting events with scope) and compared the obtained flux with previous mea- 
L > 0.4, which is indicated by the arrow in Fig 3. The re- surements (Tanimori et al, 1998b; Aharonian et al, 2000). 
suiting a distribution is plotted in Fig. 4. The data points They are consistent with each other within 12% at ~ 3 TeV. 
with statistical error bars were obtained by ON-source run. On the other hand, the statistical error of the Crab measure- 
and the hatched histogram by OFF-source run, respectively. ment was 17%. In order to check cut dependence of the 
The number of excess events was 946±108(8.7c7). In order analysis, we also changed the cut value of L, the cluster- 
to check the algorithm, we tried changing the background ing methods, and the threshold, and obtained the fluxes each 
sample, i.e., from OFF- to ON-source events to make proton time. From the above procedures, we estimated a systematic 
PDFs. The number of excess eventschanged to 938±108(8. 7(7), uncertainty of 15.4%. The stable region of L in the accep- 
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tance calculation is shown in Fig 3 by the horizontal arrow. 
The acceptance is considered to be stable within 6.5% in this 
region with respect to the L cut. 

3.1 Energy Spectrum 

For the energy scale ambiguities, we considered the error on 
mirror reflectivity, PSFs, the effects of Mie scatterings, and 
the ambiguity on the single-photon pulse height. The esti- 
mated value was 20%. 

3.2 Angular Resolution 

The pointing error of the telescope system can be checked 
by looking at bright stars in various observation periods. To 
date, we have only had bright stars beyond one-degree from 
the centre of the field of view where there were significant 
edge effects and aberrations. For now we have adopted a 
systematic error of 0.1 degree. 

4 Discussion 

An improvement in the S/N-ratio was clearly demonstrated 
so far. Generally speaking, the gamma-ray signal increased 
with the background remaining the same, when compared 
with the traditional square cut. 

As for the asymmetry parameter, only a positive or neg- 
ative cut can be applied in the square cut analysis, greatly 
reducing the number of events accepted. However, a likeli- 
hood analysis can save this situation. Almost aU of the pa- 
rameters can be put into the PDFs, even if there are small 
differences between gamma-rays and protons. It is neces- 
sary, however, to be careful of the dependencies between the 
various parameters. 

The most important thing is that this kind of analysis can 
reduce any "human bias", because the number of cut param- 
eters is very small. The systematic error, therefore, should 
be smaller than the square cut analysis. The small number 
of cut parameters makes it easy to estimate the systematic 
errors. 

Usually, the cut value of L should be around ~ 0.5; allow- 
ing automated analysis for any situations such as different 
elevation angles, etc. For the analysis of stereoscopic ob- 
servations (Enomoto et al., 2001), the product of L of many 
telescopes can be used as a single-cut parameter. In this case 
it is necessary to tune the Monte-Carlo simulations as accu- 
rately as possible. 



ysis. Also we confirmed the systematic error of this analysis 
is sufficiently small. 

Acknowledgements. This work was supported by the Center of Ex- 
cellence, and a Grant-in-Aid for Scientific Research by the Japan 
Ministry of Education, Science, Sports and Culture. Also it was 
supported by Australian Research Council. 

References 

Aharonian, F. A., et al., ApJ, 539, 314-324, 2000. 
Ellison, D. C, Berezhko, E. G., and Baring, M. G., ApJ, 540, 292, 
2000. 

Enomoto, R., et al. Astropart. Phys., in press. 

Kawachi, A. et al., Astropart. Phys. 14, 261, 2001. 

Hara, S., et al., these proceedings, 2001. 

Hillas, A. M., Proc. 19th ICRC, 3, 445, 1985. 

Hillas, A.M., et al., ApJ, 503, 744, 1998. 

Kifune, T et al., ApJ, 438, L91, 1995. 

Koyama, K., et al., PASJ 49, L7, 1997. 

Mastichiadis, A., A&A, 305, L53, 1996. 

Mastichiadis, A. and de Jager, O. C, A&A, 311, L5, 1996. 

Mori, M., et al., these proceedings, 2001. 

Muraishi, H. et al., A&A, 354, L57, 2000. 

Naito, T., et al. Astron. Nachr., 320, 206, 1999. 

Pfeffermann, E., and Aschenbach, B. 1996, in Roentgenstrahlungn 

from the Universe, ed. H. H. Zimmermann, J. Trumper, and H. 

Yorke (MPE Rep. 263; Garching: MPE), 267. 
Pohl, M., A&A., 307, 507, 1996. 
Slane, R, et al., ApJ, 525, 357, 1999. 
Tanimori, T. et al., ApJ, 497, L25 and Plate L2, 1998. 
Tanimori, T. et al., ApJ, 492, L33, 1998. 
Tanimori, T., et al., these proceedings, 2001. 

Yoshida, T. and Yanagita, S., Proc. 2nd INTEGRAL Workshop 'The 

Transparent Universe', ESA SP-382, 85, 1997. 
Yanagita, S. and Nomoto, K., Proc. 3rd INTEGRAL Workshop 

'The Extreme Universe' , Astrophys. Lett, and Comm., Vol. 38, 

461-464, 1999. 
Weekes, T. C, et al. ApJ, 342, 379, 1989. 



5 Conclusion 

We have measured gamma-rays from RXJ1713— 39 in the 
energy range between 400 GeV and 5 TeV. The statistical 
significance of the measurement is greater than 8 a. We have 
used the new likelihood method of analysis. We obtained a 
better signal-to-noise ratio compared with the standard anal- 



